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ABSTRACT 


A  multi-step  approach  to  rectify  radar  images 
geometrically  and  radiometrically  is  presented.  The 
method  uses  digital  elevation  models  (DEMs)  and  an 
image  simulator  to  establish  the  relation  between  radar 
image  coordinates  and  ortho  image  locations.  The 
simulation  serves  to  create  a  synthetic  gray  value 
image  similiar  in  its  geometry  to  the  real  radar  image, 
associating  with  each  image  pixel  a  backscatter  value 
dependent  on  the  incidence  angle  and  a  DEM  address. 
Correlation  between  simulated  and  real  image  results  in 
a  relation  between  radar  backscatter  values  and  DEM 
addresses  and  produces  an  ortho-image  in  an  additional 
step. 

Stereo-derived  and  map-derived  elevation  models 
were  used  to  demonstrate  the  rectification  process  for 
air-  and  spaceborne  radar  images. 


1.0  INTRODUCTION 


Numerous  remote-sensing  applications  are  based  on 
precisely  relating  the  gray  values  of  an  aircraft  or 
satellite  image  to  the  corresponding  phenomenon  on  the 
ground,  to  gray  values  of  other  images,  or  to  symbolic 
descriptions  of  an  object  in  the  form  of  maps. 

One  concept  to  solve  this  correspondence  problem  is  by 
physically  generating  a  derived  image  that  is  geometrically 
corrected  and  therefore  has  the  geometric  properties  of  a 
map.  The  result  can  be  denoted  by  "ortho-images"  or  by 
"geo-coded  imagery".  The  first  expression  results  from 
standard  photogrammetr ic  terminology  where  since  about  1955 
the  ortho-photo  has  become  a  routine  product  (Manual  of 
Phctogrammetry ,  1994).  "deo-coding"  is  a  more  recent 
terminology  that  was  independently  developed  in  the 
remote-sensing  and  image  processing  fields  (Guertin,  Shaw, 
1981  )  . 


The  generation  of  digital  geo-coded  ortho-images 
requires  that  given  raw  input  images  are  resampled  in  a 
digital  image  processing  system.  Resampling  must  result  in 
a  new  output  image  so  that  the  image  gray  value  of  each 
imaged  object  is  presented  in  its  proper  position. 
Essentially  therefore  one  has  to  shift  input  image  positions 
x,  y  into  output  positions  X,  Y  such  that 

-  effects  of  the  sensor  position  and  attitude  during 
image  acquisition; 

-  effects  of  perturbations  of  the  internal  sensor 
operations  ; 

-  effects  of  the  imaging  characteristics  or  type  of 
geometric  projection; 

-  effects  of  the  object  shape 

are  all  eliminated  from  the  image  geometry.  The  resulting 
output  image  has  a  specified  geometry  that  is  independent  of 
the  specifics  of  sensor  and  image  acquisition. 

This  report  describes  one  procedure  to  take  a  digital 
radar  image,  recordings  of  the  position  and  attitude  of  the 
sensor  platform,  and  a  description  of  terrain  shape  in  the 
form  of  a  digital  elevation  model  (DEM),  to  obtain  a  digital 
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ortho-radar  image.  The  approach  is  based  on  a  pixelwise 
re-arrangement  of  the  input  data;  therefore  one  can  refer 
to  the  procedure  as  "differential  rectification". 

As  a  by-product  of  geometric  rectification  one  can  also 
manipulate  the  image  gray  values  themselves  in  a  step  of 
"radiometric  rectification",  thereby  eliminating  image  gray 
value  effects  of  topographic  terrain  slopes. 

The  report  will  first  review  the  current  understanding 
of  radar  rectification  and  single  image  radar gramme t ry .  It 
will  then  discuss  the  individual  steps  of  the  rectification 
procedure  developed  in  the  current  context,  and  finally 
demonstrate  results  using  aircraft  and  satellite  images, 
both  with  map-derived  and  with  stereo-radar  terrain 
elevation  data. 
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REVIEW  OF  SINGLE  IMAGE  RADARGR AMMETR Y  FOR  RECTIFICATION 


2 . 1  Image  Coordinates 

Radargrammetry  requires  first  of  all  that  the  radar 
measurements,  time  t  and  slant  range  r,  are  properly 
extracted  from  the  radar  recording.  This  implies  some 
transformation  of  an  image  coordinate  pair,  x,  y  or  pixel 
line  and  sample  numbers,  i,j,  into  r  and  t.  The  basic 
concept  is  illustrated  in  Figure  2.1. 


We  measure  for  each  image  point  P  its  xioYl< 
comparator  coordinates  or,  if  we  deal  with  a  digital  image, 
its  corresponding  line  and  sample  numbers.  These  need  to  be 
converted  to  the  time,  t,  when  object  P  was  imaged  (sensor 
events  time),  and  range,  r,  to  P.  Numerous  methods  exist  to 
relate  xk>Yk  to  r » t .  CLenerally  we  may  therefore  just 
define  an  operator  F : 


r  =  VXk’  yk} 
C  =  FtUk’  V 


('2.1) 


FIGURE  2.1:  Measuring  coordinates  in  a  side¬ 
looking  radar  image;  xk>Yk  are  compa¬ 
rator  coordinates,  r,t  are  range  and 
time  coordinates. 


In  certain  types  of  radar  images,  particularly  those 
designed  for  good  geometry,  one  may  find  auxiliary  reference 
marks,  such  as  those  of  Figure  2.1,  denoted  by  R  ■]  ,  R£ 
etc.  Depending  on  the  type  of  marks  we  can  define  the  form 
of  Equ.  (2.1)  and  its  parameters.  As  an  example,  a  single 
row  of  so-called  time-marks  may  enable  one  to  define  a 
straight  line  to  be  denoted  as  "image  time  coordinate"  axis, 
x . 


Numerous  methods  could  be  defined  to  actually  implement 
operators  F p ,  F t .  Factors  to  be  considered  are  the  type 
of  geometric  referencing,  the  type  of  radar  record  in  a 
slant  or  ground  range  presentation ,  the  individual  digital 
or  electronic  pre-processing  of  the  image  data  etc. 

In  most  cases,  radargrammetry  begins  with  the 
SLR-image,  not  with  the  raw  electronic  signals  acquired  by 
the  antenna. 

However,  a  trend  in  digital  radar  systems  is  developing 
to  be  beyond  the  radar  record  and  to  look  into  raw, 
unprocessed  signals.  First  indications  are  visible  in 
studies  by  Kratky  (1979),  Brown  et  al.  (1981)  and  Curlander 
et  al.  (1981).  These  new  approaches  rely  on  concepts  such 
as  "sensor  events  time",  i.e.  the  exact  time  when  a  raw 
pulse  was  transmitted  from  and  received  by  the  radar 
antenna,  and  how  the  raw  radar  echo  is  converted  into  image 
points  . 

Data  from  the  SEASAT-SAR  experiment  can  serve  to 
illustrate  the  problem  of  extracting  t,  r  from  the  image. 

For  a  long  time  the  digitally  correlated  SEASAT-SAR 
images  did  not  show  any  time  or  range  reference  marks.  The 
user  was  not  provided  with  means  to  enable  him  to  relate,  x, 
y  to  t,  r  coordinates.  However,  Curlander  ( 1 9 8 1  )  recently 
established  a  relationship  that  exists  between,  on  one  hand, 
the  sensor  events  time  t  and  slant  range  r  and  on  the  other 
hand,  the  image  pixel  line  (azimuth)  number  i  and  sample 
(range)  number  j.  This  relationship  is  for  digitally 
correlated  SEASAT-SAR  rather  complex.  This  is  illustrated 
by  the  formulas  that  were  presented  by  Curlander  (  1981),  who 
went  back  to  the  raw  signal  data  to  define  r,  t  for  each 
pixel . 
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2 • 2  Projection  Equation 

Numerous  authors  have  formulated  rigorous  radar 
projection  equations,  similar  to  the  so-called  rigorous 
photogrammetric  equations  for  photography.  The  intent  of 
these  formulations  is  to  model  the  typical  range  projection. 
Object  points  are  projected  along  concentric  circles  onto  a 
projection  plane  (Figure  2.2).  Authors  who  have  contributed 
here  are  Rinner  (1948),  Konecny  and  Derenyi  (1966), 
Rosenfield  (1968),  Akowetzki  (1968),  Derenyi  (1970),  Gracie 
et  al.  (1970),  Leberl  (1970),  Hockeborn  (1971),  Norvelle 
(1971),  Geier  (1972),  Greve  and  Cooney  (1974),  DBA-Systems 
(1974)  and  many  authors  since.  A  review  was  compiled  by 
Leberl  (1976b). 
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FIGURE  2.2:  Range  projection.  Each  image  point 

defines  a  circle  as  the  locus  of  the 
corresponding  object  point. 

For  such  a  formulation  we  need  to  define  the  sensor 
position,  s(t)  =  (sx(t),  s  y ( t )  ,  sz(t))T,  and  antenna 
attitude  with  J(t),  cj  ( t )  ,  <(t).  In  synthetic  aperture 

radar,  the  synthetic  antenna's  attitude  is  merely  the  first 
derivative  of  position,  namely  the  velocity  vector  s(t)  = 
(sx(t),  s„(t),  sz(t)).  Both  s(t)  and  s.(t)  are 

functions  of  time  t.  The  formulation  also  requires  a 
definition  of  object-  and  sensor  coordinate  systems 

according  to  Figure  2.3  with  unit  vectors  ,  e_2,  e^ 

and  u,  v,  k,  respectively. 

We  then  find  for  an  object  point  P  its  position  vector 

£ : 

£  =  Pxi,  +  Pyl2  +  Pz-3  (2-2) 

and  in  the  sensor  system 

r=uu+vv+ww  (2.3) 

—  p—  p—  p—  J  ’ 

< 
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FIGURE  2.3:  Definition  of  object  space  anc 
sensor  coordinate  systems. 


Iheref ore , 


p  =  s  +  r 


t  2  .  4  ) 


It  is  straight  forward  to  see  from  Figure  2.4  that 

u  =  r  sin: 

P 

1  ■>  |  y 

vp  =  r  sin'fi  -  sin-:)  -  (2.5) 

v  =  -  r  cosfi 
P 

Auxiliary  angle  t  is  a  system  constant  called  squint-angle 
by  most  authors  (Rosenfield,  1968;  Hockeborn,  1971)  and 
defines  the  conical  shape  of  the  radar  pulse. 


Introducing  a  new  vector  £  =  (up,  vp,  wp)  in  the 
sensor  coordinate  system  and  a  rotation  matrix  A  to  rotate 
the  sensor  system  £,  v,  w  into  the  object  system  e ;  ,  e  o  , 


we  find  in  matrix  notation 


-3  » 


FIGURE  2.4:  Relating  squint  angle  ~ ,  slant  range 
r  and  elevation  angle  .1  to  sensor 
coordinates  u,  v,  w. 


where 


p  =  _s  +  A  •  p 


!  a 

( t ) 

u  (  t  ) 

u 

It)! 

1  x 

y 

2 

| 

A  =  '  v 

u) 

v  (t) 

V 

(t): 

i 

—  X 

y 

7 

'  w 

( c ) 

w  (  t ) 

W 

( c  )  1 

I  x 

L 

V 

2 

J 
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Rotation  matrix  A  can  contain  the  classical  orientation 
angles  i  ,  oj  and  <  if  we  deal  with  a  real  antenna.  With 
synthetic  antennas,  we  obtain: 


(2.8) 


u  = 

(u  ,u  ,u  ) 

=  s  / 

s  : 

X  V  z 

V 

'v  ,v  ,v  ) 

=  (_s 

x  s )  /  ;  s  x  s 

— 

x  y  z 

w  = 

(w  ,w  ,w  ) 

=  (u 

x  v)  /  |  u  X  V 

— 

x  y  z 

The  unit  vector  u  extends  along  the  synthetic  antenna,  thus 
in  the  direction  of  the  velocity  vector  s .  The  unit  vector 
v  should  be  normal  to  u  and  to  the  nadir  direction  of  the 
antenna.  For  satellite  sensor,  the  nadir  direction  is 
easily  seen  to  be  the  direction  of  the  sensor  position 
vector  in  a  planetocentric  system.  With  aircraft  radar,  one 
has  to  specify  here  explicitly  that  a  geocentric  system  is 
assumed  also  for  the  sensor  position  vector  £  and  the  object 
coordinate  system  defined  by  e-j,  e2  ,  e^. 

Equations  (2.5)  and  (2.6)  relate  object  and  image, 
where  r  and  time  t  derive  from  the  image  coordinates.  We 
now  see  that  -  ''ontrary  to  classical  photogrammetry  -  we  do 
not  algebraically  relate  image  and  object  coordinates,  but 
the  image  only  defines  time,  which  in  turn  defines  the 
parameters  s  and  A  of  Eq.  (2.6). 

If  the  imaging  arrangement,  and  thus  _s  (  t )  ,  A(t),  are 

known  (measurement  of  navigation  data),  then  we  have  for 
each  point  3  equations  ( Eq .  2.5)  with  9  unknowns  i.e.  2  , 

Px,  Py  and  pz.  It  may  be  preferable  to  work  with  2 
equations  and  3  unknowns.  This  is  achieved  by  elimination 
of  the  unknown  elevation  angle  12  (angle  between  nadir  and 
line  of  sight),  thus  giving  the  following  two  equations: 

_p  -  s_!  =  r  (2.9) 

u  •  (jp  -  s)  =  sin:  .  u  .  ,  £  -  £  =  sinx-r  (2.10) 

where  we  now  have  a  simple  geometric  description  of  the 
radar  projection  case.  The  geometric  locus  of  a  point  is 
the  intersection  of  range-sphere  (Eq.  2.9)  with  a  cone  (Eq. 
2.10).  This  is  the  circular  projection  line  mentioned 
before. 

In  Eq .  (2.10)  the  sensor  coordinate  axes  are  all 
described  by  the  geometric  position  and  velocity  vector  of 
the  antenna.  This  applies  in  particular  for  SAR,  since  the 
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antenna  attitude  is  defined  by  the  velocity  vector,  (the 
"synthetic"  antenna  extends  in  the  direction  of  the  physical 
antenna's  motion).  For  real  aperture  radar,  the  real 
antenna's  attitude  itself  is  to  be  used  to  define  u ,  v  and 
w.  Therefore  angles  b  and  <  must  be  entered  in  Eq.  (2.7); 

a j  is  insignificant  and  can  be  set  at  zero.  This  is  obvious 
from  Eqs.  (2.9)  and  (2.10)  which  do  not  contain  an  angle  u> 
since  the  components  of  vector  u  in  the  rotation  matrix  A, 
namely  ux,  uy  and  u2  do  not  contain  w.  The  same  fact 
is  also  obvious  from  geometry:  the  projection  circle  is  not 
affected  by  a  sensor  roll,  w. 

In  aircraft  radar,  r  normally  is  zero,  so  the  cone  of 
Eq.  (2.10)  degenerates  to  a  scanning  plane: 

Li-(_p-_s)=0  (2.11) 

The  range  sphere  has  an  obvious  physical  explanation:  the 

echo-time  measurement  of  radar  clearly  associates  with  each 
image  point  a  slant  range  r,  and  therefore  a  sphere,  as  one 
geometric  locus  for  the  corresponding  object. 


2 . 3  Inner  Orientation 

The  parameters  of  the  relationship  between  image 
coordinates  x,  y  and  the  object  space  entities  r,  t  are 
denoted  as  parameters  of  interior  orientation.  Equation 
(2.1)  was  a  mathematical  representation  of  this  interior 
orientation.  The  concept  derives  from  photogrammetric 
cameras,  where  each  image  point  defines  a  projection  ray 
attached  to  the  camera.  The  exterior  orientation  relates 
the  projection  ray  to  the  object  point. 

However,  Eq .  (2.1)  represents  merely  a  model  of 
reality.  Errors  of  this  model  must  be  anticipated  as  shown 
in  Figure  2.5.  Of  course,  any  error  model  has  unknown 
coefficients  and  the  question  is  how  to  obtain  estimates  of 
these  coefficients.  As  will  be  shown  later,  these  can 
derive  from  known  object  point  coordinates.  Error  models 
generally  imply  a  polynomial  form  of  the  image  coordinate 
deformation  Ar ,  At.  Gracie  et  al.  (  1  970  ),  Leberl  (1971), 
DBA-Systems  (1979)  and  others  have  applied  the  following 
polynomials : 


information  to  relate  image  and  object  space.  This  has  been 
attempted  by  many  authors.  The  efforts  can  be  classified 
into: 

a)  interpolative  models  employing  so-called  "rubber¬ 
sheeting"  or  warping  algorithms  to  make  the  image 
fit  to  ground  control; 

b)  parametric  models  employing  projection  equations 
and  computing  unknown  parameters  from  control 
point  information;  and 

c)  hybrid  combination  of  a)  and  b),  where  the  pro¬ 
jection  equations  are  used  to  relate  image  and 
object,  but  where  an  interpolative  technique  is 
applied  to  represent  unknown  deviations  of 
reality  from  the  assumed  model. 


2.4.1  I  n  t,  e  r  do  la  t  i  ve  Methods 

A  first  approach  to  the  study  of  image  geometry  is  to 
consider  the  image  itself  as  a  map  substitute  and  to  analyse 
its  differences  from  a  map  of  known  geometric  properties. 
Differences  ^Xj_,  Ay^  between  image  Xj_,  y and  map 
xi>  *i  coordinates  can  be  measured  at  control  points  i. 

Some  interpolation  function  is  then  used  to  predict  Axp, 

A y p  deformations  of  the  image  at  arbitrary  image  locations 
p.  The  same  methods  can  be  employed  as  those  used  in  other 
photogrammetr ic  interpolation  tasks  (Leberl,  1975).  To  give 
but  two  examples,  let  us  consider  Figure  2.6.  The  points 
?i,  i=1,...,n  are  known  with  their  image  coordinates  Xj_, 
y  i ,  and  map  coordinates  X  ^ ,  Yj_. 
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(a)  Polynomials 


We  may  now,  as  one  example,  attempt  a  scale  fit: 


rx- 

al1  a  1  2 

X  . 

1 

Y  . 

l 

a 1 2  ai1 

_yi_ 

FIGURE  2.6:  I n t e r po la t i v e  method  of  relating  image 
and  map  points;  a  preliminary  trans¬ 
formation  leads  to  residuals  AX,  AY. 
These  define  a  correction  surface  that 
can  be  described  by  isolines,  separa¬ 
tely  for  AX  and  for  AY. 


Since  this  is  solved  with  more  control  points  than 


needed 


one 


obtains  residuals  AX,  AY  or,  if  Eq .  (2.13)  is 


find  thus  first  Xp, 


Ay. 

yp 

For  a  new  point,  xp,yp 

we 

and  we  obtain  the  final 

map 

AXp, 

AYp: 

X '  + 

AX 

p 

P 

(2 

Y  '  + 

AY 

P 

P 

(2.14) 


where  AXp,  AYp  result  from: 


AX 
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a  +  a 
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X  +  a 
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Y  +  a 
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,X  Y  +  a  Y 
3  p  p  4 
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(2.15) 


pr-T  w  T  ^ 
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The  coefficients  a0,...,an,  bg,...,bn  are  computed 
from  the  known  control  point  values  AXj_,  Ay^_. 

Various  authors  have  evaluated  different  types  of 
polynomials  for  radar  image  correction  (Leberl,  1970; 
Derenyi,  1974a;  Hirsch  and  van  Kuilenburg,  1976;  Shakine 
and  LeToan,  1978).  A  common  conclusion  is  difficult  to 
draw,  since  different  data  from  different  sensors  were 
employed.  However,  one  may  summarize  that  the  use  of 
polynomials  is  generally  not  advisable  due  to  their 
predictably  unfavourable  behavior  in  the  case  of  irregular 
control  point  distribution,  and  in  the  case  of 
extrapolation.  Polynomials  are  only  recommended  if  they  are 
of  low  order  (1st  or  second)  and  if  there  is  an  even  control 
point  distribution. 


(b)  Moving  average 

The  polynomials  of  Eq.  (2.15)  are  a  global,  closed 
representation  of  image  deformation.  There  are  also 
point-wise  methods  of  describing  this  deformation.  We  start 
from  the  known  set  of  A  X  ,  AYj_-values  in  control  points, 
and  need  to  compute  the  A  X  p ,  A  Yp  in  a  new  point  Xp, 

Yp.  The  moving  average  technique  assembles  for  each  point 
P  a  set  of  known  control  points  according  to  some  criterion. 
This  may  for  example  be  the  n  closest  control  points  (Figure 
2.7),  or  it  may  be  m  control  points  in  each  of  four 
quadrants  (Figure  2.8). 
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'IGURE  2.7:  Moving  average  interpolation  with 

square  windows  to  collect  n  reference 
points. 


FIGURE  2.8:  Selection  of  control  points  in 
quadrants,  in  this  example  with 
2  points  per  quadrant. 

With  these  selected  control  point  sets  a  polynomial  of 
low  order  is  now  computed  (1st  and  2nd  order)  whereby  each 
control  point  obtains  a  weight  as  a  function  of  its  distance 
from  the  new  point.  If  all  coordinates  have  been  reduced  to 
an  origin  in  the  new  point,  then  only  the  first  coefficients 
&q,  bg  of  this  polynomial  need  to  be  computed: 
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AX  =  a  +•  a  X '  +  .  .  . 
P  o  Ip 

AY  =  b  +  b,Y'  +  ... 
P  o  Ip 


since  X',  Yp  =  0. 

This  computation  and  control  point  selection  is 
repeated  for  each  new  point.  In  the  case  of  a  first  order 
polynomial  one  also  has  the  expression  "moving  plane" 
interpolation.  A  zero-order  polynomial  leads  to  a  simple 
weighted  average  i n t er pola t ion . 


2.4.2  Parametric  Models 


The  exterior  orientation  of  the  sensor  is  a  function  of 
time  and  therefore  continuously  changing.  As  a  result  each 
image-x-coordinate  relates  to  a  different  exterior 
orientation.  It  is  obvious  that  there  cannot  be  rigorous 

reconstruction  of  the  exterior  orientation  from  control 

points  for  all  times  t:  the  number  of  unknown  exterior 

orientation  elements  is  infinite.  Even  if  we  had  several 

control  points  imaged  at  one  given  time  t  and  thus  along  one 
line  in  the  image-y-direction,  complete  computation  of  the 
exterior  orientation  is  not  possible. 

To  illustrate  the  problem  let  us  assume  that  the 
exterior  orientation  is  nearly  perfect  and  that  the  flight 
only  slightly  deviates  from  a  straight  line  along  the  object 
e  i -coord  ina  te  axis.  Also  <P  ,  ui  ,  <  ,  the  elements  of  the 
sensor  attitude,  are  nearly  zero. 

These  small  deviations  dj,  da,  dx,  dsx, 
dsz  of  the  sensor  from  a  predictable  exterior 
can  be  related  to  image  coordinate  errors,  dr, 
dy.  We  start  from  Eqs.  (2.5)  and  (2.6): 

£  =  A,-(£-_s)  (2.16) 

Differentiation  of  Eq.  (2.16)  leads  to: 


dSy  and 
orientation 
d  t ,  or  d x  , 
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where  d£  =  (dsXf  dsy,  dsz).  Since  we  assume  0  =  w  =  <.*o, 
we  get  A-  to  be  the  unit  matrix,  E. 

The  differential  form  d A- ^  •  (£  -  s)  is  well  known 

in  photogrammetry  (Rinner  and  Burkhardt,  1972): 
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[~  0  -d< 

d<t>i 

d  A_  '  ■ 

s )  = 

;  d  K  0  - 

d  oj  |  A 

1  (£"£.) 

(2 

| _ —  d  <t>  d  in 

oj" 

we  get 

with  A~  1 

(£  -  s) 

= 

(  u  p  >  v  p  > 

!~du  i 

i  P 

!  -d<  • 

v  +  d 0 • w  -ds 

P  P  x 

t 

;  d  V 

1  Pi 

! 

d  <  •  u 

-du • w  -d  s 

P  p  y 

(2 

(2.18) 


LdWpj  L~~'  'p  "zJ 


•  d  0  •  u  +doj.v  -ds 


Errors  of  image  x , y-coordinates  result  from  duD,  dvD  and 


dWp  as  f  ollows  : 


d  x  =  du 


But  y=r= ( vp2+Wp2 ) 1 /2 #  Therefore: 


dy  =  dr  =  —  (v  • dv  +w  • dw  ) 
r  p  p  p  p 


=  sinft-dv  -cosft’dw 
P  P 


Let  us  now  simplify  further  with  squint  angle  t=0. 
up=0,  Vp=r*sinft  and  Wp=-r*cosfi.  Thus 


Then 
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dv  =  si  nil  (-duj  -w  -ds  )  -  cosiHduj-v  -ds  ) 

P  Y  p  z 


dv  =  ♦  du  sinQ1 

1  r • cosQ  - 

s  inQ 

d  s 

V 

-  c  o  s  .1  d  -  ■ 

1  r  *  s  i  n  i*  + 

COS-”'. 

d  s 

(2.20) 
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dv  =  -  sin:',  ds  +  c  o  s  2  ds 
y  z 

dx  =  -  r  ‘  sin.'.  d<  -  r- cos-1  di  -  ds 

x 

With  flying  height  H  =  r*cos:l,  or  with  object  coordinates 
Py=r*sinft,  this  can  be  rewritten: 

dx  =  -  ds  +  p  •  d  <  -  H  •  do 
x  v 


y  y  z 

We  see  again  that  we  have  only  five  elements  of  exterior 
orientation.  The  sixth  element,  d  ,  has  no  effect  on  image 
geometry . 

For  each  image-x-coordinate,  a  complete  set  of  five 
orientation  elements  must  be  solved.  This  is  not  possible 
with  a  co-linear  arrangement  of  control  points,  since 
clearly  the  effects  of  dsx  and  d  d>  cannot  be  separated.  We 
can  solve  either  (dsx,  dsy,  ds2,  d<)  or  (do,  dsy, 
ds2,  d<).  The  fifth  element  must  be  known.  The  required 
number  of  control  points  to  compute  the  four  elements  is 
two . 

It  is  unrealistic,  however,  to  assume  that  for  each 
image  x-line  one  can  have  two  control  points.  Also  this  is 
not  necessary  because  of  correlation  that  exists  between  the 
exterior  orientation  values  at  two  different  times  t-|  and 
t2  . 


Authors  have  in  the  past  proposed  various  mathematical 
models  for  the  exterior  orientation  as  a  function  of  time, 
e.g.  DBA-Systems  (197*0  employed  polynomials  such  that: 

7 

dp(t)  =  a  +  a  t  +  a  t  “  +  ... 

o  1  2 

=  f  +f.t+f,t2+... 


d  s  ( t  ) 


(2.22) 


Instead  of  solving  d  ?  ,  d<,  d  s  x ,  dsy  and  dsz  one  solves 
for  unknown  polynomials  coefficients  aQ, 

a1,...  o*  !»•••  Eq.  (2.21)  is  thus  rewritten: 


c  +c.t+c.t 
o  1  2 


+  p  b  +  b  .  t  +  b  -  t  +  .  . 
v  o  1  2 
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(2.23) 


—  ‘id  +d,t+d-t“+... 
r  o  1  2 


H.  •  2 

—  ■ e  +e.t+e,t  +...; 
r  o  1  2 


Not  all  coefficients  can  be  determined  separately.  The 
analysis  of  this  approach  is  beyond  the  scope  of  this  paper. 
It  may,  however,  be  relevant  to  point  out  that  instead  of 
single  high  order  polynomials  it  is  preferable  to  use 
piecewise  polynomials  of  low  order.  Baker  and  Mikhail 
(1974)  have  studied  these  approaches,  however,  for  aircraft 
scanner  imagery,  not  for  radar. 

Dowideit  (1977)  studied  the  method  of  modelling  the 
exterior  orientation  parameters  by  Fourier  series,  correctly 
arguing  that  the  elements  of  the  exterior  orientation  have  a 
periodical  behavior  with  several  dominating  frequencies 
(Leberl,  1972a).  Again,  one  has  to  solve  for  unknown 
coefficients  of  the  Fourier  series.  The  main  problem  is 
instability  of  the  resulting  equation  system  due  to 
correlation  among  coefficients. 

A  somewhat  different  approach  has  been  proposed  for 
dynamic  scanner  images  by  Ebner  and  Hoessler  (1978).  The 
elements  of  exterior  orientation  are  considered  to  form 
time-series  appropriately  modelled  by  Gauss-Markof f-chains 
of  order  m: 


MCn+l) 


a  .  s  ( t  .  )  +  e 

l  n-m+ i  n+ 


(2.29) 


Thus  the  exterior  orientation  is  described  at  discrete  times 
b i , t2 , . . • 5  at  time  tn,  the  orientation  parameter  is 


considered  to  be  a  linear  extrapolation  from  m  previous 
times  plus  a  random  addition,  e  .  Ebner  and  Hoessler  could 
show  that  a  numerically  stable  and  comparatively  manageable 
solution  can  be  formulated.  A  detailed  description  is 
beyond  the  scope  of  this  paper.  An  application  to  radar 
must  still  be  investigated. 


2.4.3  Hybrid  Methods 

A  combination  of  parametric  and  interpolative  methods 
is  feasible  and  common  in  photogrammetry.  A  mathematically 
rigorous  model  of  the  imaging  process,  such  as  used  in  the 
parametric  case,  may  not  be  able  to  describe  the  entire 
chain  of  physical  events  leading  to  the  image.  Unknown 
elements  or  effects  must  be  considered.  This  can  be  done  by 
introducing  Ax,  Ay-corrections  to  the  image  coordinates; 
these  corrections  can  oe  represented  by  polynomials  or  other 
interpolative  methods. 

In  photogrammetry  this  approach  is  denoted  by 

"self-calibration",  the  concept  is  related  to  the  method  of 
so-called  "collocation"  (Moritz,  1973)*  An  elementary 
implementation  is  by  solving  first  a  parametric  formulation. 
This  will  be  done  in  a  least-squares  method  that  leaves 
residuals  at  the  observed  image  points.  The  residuals  are 
then  subject  to  an  interpolative  computation.  This  is 
meaningful  only  when  the  residuals  form  a  systematic 
pattern.  A  new  point  P ( x p  ,  y 0 )  is  thus  converted  to 
object  coordinates  by  first  adding  Axp,  Ayp  to  xp, 
yp ,  and  by  then  applying  the  radar  equations. 


2 . 5  Projection  Of  A  Ground  Control  PoiHi,  Into  A  Radar  Image 

In  the  event  that  the  exterior  orientation  is  measured, 
e.g.  with  the  help  of  inertial  navigation,  one  must  expect 
errors  of  measurement  to  exist.  These  require  correction 
with  the  help  of  ground  control  points. 
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Figure  2.9  illustrates  the  solution  to  the  computation 
the  image  coordinates  x,  y  for  a  given  object  point 
,  ?v>  p 7 )  as  it  was  published  by  Leberl,  Fuchs  and 

(  1  9  8  1  )  .  The  exterior  orientation  is  known  in  the  form 
position  vectors  Sj_ }  i=1,...,m  at  times  tj_, 

.  .  .  ,  m .  Also  given  is  thus  the  velocity  vector  _s  ^  = 

-  s j_ ,  or  the  equivalent  attitude  vector  describing 
,  <Tt ) . 


FIGURE  2.9:  Computation  of  flight  positions  ^  p 
at  which  ground  point  P  was  imaged. 
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From  figure  2.9  one  finds: 

(s.  -  p  )  •  s  .  =  s.  -  p  •  s.  cosa. 

—  i  —  —  l  —l  —  —  L  1 

where  a  j.  is  'he  angle  between  the  linear  orbit  element  i 
and  the  vector  from  Sj_  to  .  It  relates  to  entity 

eg  : 

cosa  =  e  /  s.  -  p 
g  -i-  ~ 


where  eg  results  from: 


(.  s  .  -  p  )  •  s  . 

—  l  —  —i 


s  .  •  e  ,  e 

-i  g  i 
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With  this  we  get  orbit  position  s  p : 


s  =  s  .  +  s  .  •  e  /  s  . 

-p  -l  -i  g  -l 


(2.25) 


Position  Sp  then  leads  to  slant  range  rp  and  time  tp 
o  f  imaging : 

r  =  (p  -  s  )  (2.26) 
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(2.27) 


Actual 
through  an 


image  coordinates  xp,  yp  are 
inversion  of  the  equations  of  inner 


obtained 
orienta  tion : 


x.+(t  -t  ) (x.  -x. )/(t  -t . ) 

i  P  i  i+i  i  l+l  i 


v  =  r  /  f  -  c 
P 

Here,  x^,  xi+i  are  image  x-coordi nates  of  time 
references  ti?  ti+1. 

These  computed  image  coordinates  can  be  compared  with 
actual  measurements  of  x,  y.  Discrepancies  are  the 
consequences  of  errors  of  inner  and  exterior  orientation: 
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X 


X 


Ax  = 

Ay  =  y  -  y 


The  Ax,  Ay -d i s c re p an c ie 5  in  control  points  can  be  used  to 
compute  correction  functions  and  then  to  compensate  for 
these  errors  in  new  points.  The  correction  functions  can  be 
of  the  types  as  described  in  section  2.4.1. 


2 . 6  Rectification 

In  orthophoto-generation  by  both  digital  image 
processing  and  photogrammetr ic  computer  assisted 
orthoprojection,  one  needs  to  define  a  sufficiently  dense, 
geometrically  correct  grid  of  anchor  points  according  to 
Figure  2.10.  These  anchor  points  are  generated  using 
rigorous  radargrammetry.  All  image  details  inside  a 
quadrangular  grid  mesh  are  transformed  with  a  simpler  and 
thus  faster  linear  method. 


FIGURE  2.10:  Anchor  points  in  the  deformed  image 
(left)  and  orthophoto  (right)  (from 
Kraus  ,  19  7  6)  . 

The  process  is  based  on  a  raster  DHM  at  orthophoto 
anchor  points.  These  terrain  points  are  resected  onto  the 
image  according  to  Eqs.  (2.25)  to  (2.27),  so  that  a 
situation  as  shown  in  Figure  2.10  results.  Based  on  these 
image  coordinates  one  of  the  standard  digital  or  analog 
p h o t og r a mm e t r i c  rectification  routines  can  be  applied. 
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Rectified  aircraft  radar  images  were  obtained  from  a 
study  by  Leberl  and  Fuchs  (1978),  and  by  Leberl,  Fuchs  and 
Ford  (1981)  with  a  detailed  error  analysis.  Its  discussion 
here  would  go  beyond  the  scope  of  this  report. 

Digital  rectification  has  been  done  so  far  in  the 
context  of  the  SE AS AT-mis s i on .  Curlander  et  al.  (  1  9 8  1  ) 
have  used  digital  SAR-images  from  SEASAT  together  with 
orbital  data,  and  Naraghi  et  al.  (1981)  matched  the 
SAR-data  with  a  digital  terrain  model. 

Instead  of  actually  creating  a  geometrically  corrected 
radar  image  one  can  think  of  converting  only  certain 
interpreted  image  details  to  a  map  projection.  If  a  DHM 
exists  with  the  image  and  the  image  is  related  to  the  ground 
height  through  a  computation  of  exterior  orientation,  then 
one  can  trace  the  image  lines  of  interest  with  a  digitizer; 
the  image  x,  y-ooordinates  are  transformed  to  X,  Y  knowing 
2.  This  concept  was  realized  by  Greve  and  Cooney  (1974) 
with  radar,  and  is  a  method  of  photogrammetry  which  has  also 
been  proposed  by  other  authors  (e.g.  Makarovic,  1  97  3  )  - 

The  idea  has  not  yet  had  responses  of  any  significance, 
either  with  metric  photography  or  with  radar  or  other  image 
types.  Interpretation  of  image  contents  apparently  relies 
on  the  use  of  the  stereo  impression  so  that  a  pure 
single-image  approach  is  not  accepted  even  with  existing 
DHMs  . 


2 . 7  PasJ;  Accuracy  Studies 

Numerous  authors  have  compared  single  radar  images  to 
maps,  thereby  employing  one  of  the  described  computing 
methods.  An  example  of  an  early  effort  was  by  Leberl 
(1972a),  as  shown  in  Figure  2.11.  Various  interpolative 
methods  were  used  to  describe  the  image  deformation. 
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FIGURE  2.11:  Example  of  a  data  set  for  an  early 
single  image  accuracy  elevation. 
British  real  aperture  radar  flight 
(  1  965  )  over  Noor d -Oos t - P o Id e r  ,  the 
Netherlands,  and  nap  at  scale 
1:50  000  for  comparison. 
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Table  I  is  a  review  of  some  existing  work.  In  view  of 
the  many  possible  mapping  methods,  of  variations  in  the 
image  types,  of  input  data  available  on  the  terrain  etc.,  it 
is  obvious  that  a  comparison  of  results  is  nearly 
impossible.  However,  one  conclusion  can  be  drawn,  the  range 
of  achievable  results  stretches  from  the  geometric  ground 
resolution  (Gracie  et  at.  1970)  to  some  large  figures. 
Clearly,  though  these  larger  figures  may  be  caused  by 
unknown  topographic  relief,  lack  of  ground  control  points  or 
poor  internal  geometry. 

Expectations  were  expressed  in  the  preparation  phase  of 
satellite  radar  that  it  could  result  in  higher  image 
geometry  performance  than  aircraft  radar  (Leberl,  1976b). 
3EA3AT-SAR  was  the  first  data  to  potentially  prove  that 
expectation.  However,  the  overall  system  engineering  for 
SEA3AT-SAR  did  not  optimize  geometric  performance.  Table  II 
presents  results  obtained  with  SEASAT-SAR.  Hardcopy  images 
on  film  either  from  optical  or  from  digital  correlation 
provided  the  same  accuracy.  Time  and  range  referencing  was 
not  of  sufficient  quality  to  fully  capitalize  on  a 
relatively  stable  orbit  geometry  and  the  resulting  good 
geometric  potential  (Leberl  et  a  1 .  1931).  Without  any 
ground  control  no  meaningful  positioning  of  ground  points 
was  feasible.  With  one  or  two  such  points  the  positioning 
error  is  +  300  to  900  m. 
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TABLE  II 

SEASAT-SAR  geometric  accuracy.  Coordinate  errors  in  km  after  radargrammetric 
measurements  and  computation,  using  the  arctic  passes  over  Banks  and  Victoria  Islands 
taken  by  SEASAT  covering  35  000  kmi  and  pass  351  over  Los  Angeles 


No.  of  control 
points 

Rigorous  radar- 
grammctry 

*  y 

Simple  scale  tit 

x  y 

Sea  arctic  images  of  Banks  Island,  35  000  sq  km 

i 

0.45 

0.49 

_ 

_ 

3  (at  both  ends) 

0.33 

0.26 

0.37 

0.37 

17 

0.17 

0.16 

— 

— 

Image  over  Los  Angeles,  7  pts  per  100  sq  km 

Optical  corr 

0.024 

0017 

_ 

— 

Digital  corr 

0.022 

0.012 

— 

— 

The  rigorous  use  of  digital  imagery  enables  one  to  more 
fully  exploit  the  data  geometrically  (Brown  et  al.  1981). 
Each  pixel  can  be  associated  with  a  precise  sensor  events 
time  and  a  slant  range  between  antenna  and  object.  As  a 
result  the  location  of  an  object  on  the  Earth's  surface  can 
be  computed  without  any  ground  control  point  to  within 
several  hundreds  of  meters.  For  the  case  of  SEASAT,  there 
is  thus  little  difference  between  the  use  of  rigorous 
radargrammetr y  and  coarse  approximations,  except  for  the 
total  lack  of  ground  control  or  the  presence  of  merely  one 
or  two  points.  For  more  appropriate  mapping  applicability  a 
satellite  radar  would  have  to  be  equipped  with  precise  time 
and  range  references  and  would  need  to  be  accompanied  by 
accurate  satellite  orbit  data  with  errors  smaller  than  the 
geometric  resolution  of  the  sensor. 
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3.0  DESCRIPTION  OF  A  SYSTEM  FOR  DIGITAL  RADAR  .IMAGE 
RECTIFICATION 


3 . 1  General 


Input  to  a  system  for  radar  image  rectification  consists  of: 

-  a  digital  radar  image, 

-  sensor  interior  orientation  parameters, 

-  sensor  position  and  attitude  values, 

-  ground  control  points, 

-  a  raster  digital  elevation  model  (DEM). 

In  the  current  context  one  has  to  use  the  input  to 
first  create  a  mathematical  relationship  between  image  and 
map  and  then  to  actually  create  a  rectified  image.  Many 
alternatives  could  be  defined.  The  system  developed  under 
this  study  is  based  on  the  assumption  that 

-  images  exist  of  mountainous  areas, 

-  identification  of  ground  control  points  is 

difficult. 

Therefore  the  system  is  based  on  radar  image  simulation  to 
support  the  generation  of  a  precise  correspondence  between 
radar  image  and  object  coordinates. 

Radar  image  simulation  with  DEMs  is  only  meaningful  in 
the  case  of  mountainous  terrain;  in  flat  areas  a  simulation 
as  used  here  would  be  redundant. 

A  simulation  system  SIMRISA  (Simulation  of  Radar  Images 
using  an  Image  Space  Algorithm)  has  been  developed  and  used 
for  radar  stereo  analyses  (Domik,  Leberl  and  Raggam,  1983). 
It  is  extended  to  serve  for  rectification  in  the  current 
effort. 
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3 . 2  Radar  Rectification  Based  On  Simulation 

The  procedure  used  for  the  rectification  of  radar 
images  of  mountainous  terrain  consist  of  the  following 
steps : 

(a)  Generation  of  a  digital  elevation  model,  either 
from  a  map  or  from  other  external  sources,  or 
from  radar  stereo  mapping; 

(b)  Verification  and  improvement,  or  computation  of, 
sensor  position  and  attitude  parameters  based  on 
identifiable  ground  control  points; 

(c)  Using  the  outputs  of  steps  (a)  and  (b),  i.e.  a 
DEM,  radar  imaging  parameters  and  a  flight  path, 
create  a  simulated  radar  image  to  resemble  the 
given  input  image; 

(d)  Correlate  the  real  and  the  simulated  radar  image, 
either  by  pointing  to  homologue  features  or  by 
automated  correlation,  thereby  creating  a  grid  of 
anchor  points  in  the  real  radar  image  to  relate 
to  the  simulated  one  and  consequently  to  the 
object  coordinate  system; 

(e)  Using  a  warping  function,  correct  the  real  image 
to  fit  over  the  simulated  image; 

(f)  Generate  a  geometrically  rectified  real  image  by 
assigning  the  real  image  gray  value  associated 
with  the  simulated  pixel  to  each  DEM  grid  cell; 

(g)  3y  substraction  of  the  real  image  gray  values 
from  the  simulated  image  one  obtains  a  radio- 
metrically  corrected  output  image. 

The  above  steps  imply  that  an  extensive  software  and  radar 
image  processing  system  is  available.  In  fact  the  system 
needs  to  consist  of  various  separate  elements.  In  the 
current  context  one  has: 

-  DEM  generation  by  the  Graz  Terrain  Model  ( GTM ) 
system  which  takes  given  map  contour  lines  to 
create  a  raster  DEM; 

-  System  for  stereo  mapping  with  radar  techni¬ 
ques  (SMART)  to  extract  contour  lines  or  raster 


DEM  data  from  stereo  radar  imagery  on  a  photo- 
grammetric  analytical  plotting  machine  (Kern 
DSR -  1  )  ; 

-  Radar  space  resection  programs  to  compute  the 
flight  parameters  from  a  given  input  image  and 
ground  control  points; 

-  Image  processing  programs  for  image  correlation 
display  and  coordinate  extraction; 

-  General  purpose  rectification  based  on  a  grid 
of  anchor  points  (RECTIF); 

-  Radar  image  simulation  (SIMRISA); 

-  Assignment  of  gray  values  to  DEM  grid  cells  based 
on  an  address  file  as  a  result  of  correlating 

the  real  with  the  simulated  image  (GEOREC). 

The  large  set  of  required  software  for  this  purpose  was 
envolved  over  the  years  in  a  sequence  of  individual  efforts. 
The  following  will  describe  some  of  the  more  relevant 
components,  in  particular  the  method  of  DEM  generation  from 
stereo  radar  images  and  the  radar  image  simulation. 

For  an  overview  of  data  and  programs  see  Table  III. 
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4.0  GENERATING  A  DEM  WITH  RADAR  STEREO  MAPPING 


4 . 1  Background 

Modern  photogrammetric  analytical  plotters  have  opened 
up  the  possibility  to  employ  other  types  of  images  than 
metric  aerial  photography  for  stereo-plotting.  Instead  of  a 
fixed  mechanical  or  optical  analogon  for  the  projection  ray 
one  has  a  flexible  computer  to  relate  image  xy-coordinates 
to  the  cartesian  X  Y  Z  object  or  model  system. 

Norvelle  (1972)  probably  was  the  first  to  actually 
program  an  analytical  plotter  for  stereo  radar.  The 
flexibility  of  earlier  analytical  plotters  was  hampered  by 
the  available  computers  and  the  use  of  low-level  programming 
languages.  Prior  to  Norvelle ’ 5  work  s tereo-radargrammetry 
was  developed  in  the  context  of  basic  studies,  e.g.  by 
Rosenfield  (1968)  and  Gracie  et  a  1  . (  1  9 7 0  )  ,  and  had  led  to  a 
rather  unique  and  specific  radar  stereo  mapping  instrument 
(Yoritomo,  1972;  Graham,  1972). 

The  current  generation  of  analytical  plotters  employs 
general  purpose  computers,  distributed  processing  and  high 
level  programming  languages.  Therefore  the  effort  required 
to  program  them  is  much  reduced  when  compared  to  the  early 
systems.  This  fact  is  combined  with  an  increased  interest 
in  radar  remote  sensing  as  documented  by  satellite  and  other 
radar  imaging  experiments.  Therefore  it  appears 
increasingly  meaningful  to  study  radar  stereo  mapping. 

In  an  effort  to  demonstrate  the  flexibility  of  current 
commercial  analytical  photogrammetric  instruments  and  in 
order  to  explore  more  fully  the  capabilities  and  limitations 
of  stereo  radargrammetry  a  procedure  and  software  system  for 
the  analytical  plotter  Kern  DSR-1  was  developed  to  set-up, 
and  plot  from,  radar  stereo  models.  This  is  a  first  and 
preliminary  report  on  this  development,  emphasizing  a  broad 
description  of  the  procedures  of  SMART  --  Stereo  M_a p p i n g 
with  Radar  Techniques.  1  ) 

The  chapter  will  discuss  the  procedure  and  demonstrate 


1)  This  development  was  done  outside  the  supported  part  of  the 
current  project  D A J A 4 5 -8  3-C -00 1 8 .  It  is  reviewed  here  due 
to  its  interest  to  the  radar  rectification  and  for  ease  of 
reference;  see  also  Raggam  and  Leberl  (1984). 


32 


it  with  a  unique  data  set.  Common  radar  image  acquisition 
is  with  parallel  flight  lines.  It  is  an  added  complexity  to 
use  intersection  flight  lines.  As  will  be  shown  the  DSR-1 
is  able  to  handle  this  type  of  data  and  produce  a  digital 
terrain  elevation  model  of  an  accuracy  commensurate  with  the 
radar  input . 


4 . 2  Overview  0_f  S^s_t£m  Modules 

The  programs  for  stereo  radar  are  organized  in  separate 
modules.  The  entire  process  of  model  set-up  follows 
radargrammetric  formulations  presented  earlier  and  reviewed 
by  Leberl  (1983).  Essentially  the  procedure  requires 
sufficient  ground  control  for  the  transformation  of  radar 
coordinates  to  the  object  space.  The  individual  images  are 
processed  first  to  obtain  approximate  values  for  the 
rigorous  radargrammetric  solution.  In  the  terminology  of 
classical  photogrammetry  one  first  solves  resections  in 
space  for  the  two  images  to  then  continue  with  a  bundle 
adjustment  using  both  images  simultaneously. 

Should  no  ground  control  be  available  then  one  has  to 
work  with  given  or  assumed  parameters  of  the  flight  path  and 
sensor  attitude.  The  preliminary  resections  would  be 
skipped . 


(  a  )  Information  ana  I  nit  i  a  1  i  z  a  t  io_n  Module:  This  task 
enables  the  input  and  manipulation  of  initial  project  and 
image  parameters.  It  serves  for  general  information  for  the 
operator,  such  as  for  example  project  name,  available  images 
for  the  project  or  image  and  orbit  data  information. 

(  b )  Control  Point  Management.  Module  :  This  module  is 
for  the  input  and  manipulation  of  ground  control  point  data. 
The  ground  coordinates  are  in  an  orthogonal  cartesian 
coordinate  system,  usually  referenced  to  an  origin  within 
the  area  to  be  mapped  and  the  Z -  a  x i s  pointing  along  the 
local  vertical  of  the  origin. 


v  c  )  Orbit  or  ?_l_ight  Data  Mans  gemejc  t  Modular  This  is 
for  input  and  manipulation  of  sensor  position  measurements. 
The  system  assumes  the  orbit  or  flight  to  be  represented 
with  time  polynomials  so  that  sensor  positions  may  be 
expressed  as  follows: 


s  =  a  + 


b  .  t  +  c  .  t 2 


where  t  is  the  time.  Approximations  for  the  o i 
coefficient  vectors  (a,  b ,  z,  ..)  will  be  determined  ir. 
orbit  data  management  module.  Input  are  the  pos i ‘ 
measurements,  output  is  an  array  of  polynomial  coefficie: 
Joint  low-order  polynomials  will  be  used  instead  of  hi. 
order  polynomials  if  the  need  exists. 


( d )  Single  Radar  Processing:  This  modui- 
represents  the  equivalent  of  a  resection  ir.  space  to  c  c  m  p  u  t  ■ 
the  orientation  parameters  for  a  single  radar  image  w i t • 
inner  and  exterior  orientation.  Plotting  with  a  s  i  n  g  1  • 
radar  image  may  be  done  on  either  a  spherical  surface  : 
chosen  radius  or  in  a  three-dimensional  X,  Y ,  Z  -  coord: rat 
system  w i t h  k town  constant  Z  . 


exterior 


(  e )  Sterejo  Radar  linage  Processing:  This 

includes  the  actual  set-up  of  a  radar  stereo  mo-; 
follows  the  use  of  the  single  radar  image  m 
determine  approximations  of  the  exterior  or 
parameters  of  the  two  overlapping  images.  After  el 
of  y-paraiiaxes  this  module  serves  to  plot  contour 
p  1  a  .nine  trie  features.  Results  are  either  directly 
on  an  xy  pen  plotter  or  are  entered  into  a  digital 
for  further  processing.  One  example  that  is 
operational  is  for  digital  terrain  elevation  models 
3  r  a  z  Terrain  Model  (3TM)  program  system  (Oswald, 
in  print). 


on  an  xy 


Radar  Stereo  Model  Set-uo 


As  previously  a  escribed  the  radar  stereo  model  set-up 


in  SMART  is  realized  in 
e  o  m  o  u  t  a  t  i  o  r.  of  the  elements 


;  a  g  e  o  a  l  r 


oliowe; 


i  two  steps  consisting  of 
:  of  inner  orientation  of 
by  an  exterior  orientation  w 


radar  bundle  adjustment. 


Inner  Orientation:  The  inner  orientation  requires  the 
establishment  of  a  relationship  between  plate  or  image  x  ,  y 
coordinates  and  the  physical  radar  measurements  of  time  t 
and  slant  range  r  (see  section  2.3).  The  system  uses  a 
so-called  'range  reference  line'  at  the  near  range  edge  of 
the  image;  it  would  correspond  to  the  start  of  the  sweep  on 
the  image  recorder.  It  should  be  defined  by  some  distinct 
.icx  marks  (fiducials).  In  the  event  that  no  such  marks 
exist  the  operator  has  to  create  artificial  marks  and 
determine  the  inner  orientation  in  a  process  of 
self -calibration  together  with  the  exterior  orientation. 


I 
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Exterior  Orientation:  Measurement  of  homologue 


orientation  points 

can  first 

of  ail 

be  as 

in  a  comparator 

by 

removing  manually 

the  x-  and 

y -par  a 11  axes . 

Secondly  after 

a 

preliminary  resec 

tion  with 

each 

image 

one  is  in  a  near 

iy 

parallax -free  stereo  model  and  observation  of  homologue 
points  is  made  more  convenient.  In  contrast  to  analog 
stereo  photogrammetry  the  analytical  plotter  does  not  need  a 
relative  orientation  task  seperate  from  absolute 
orientation.  The  parameters  of  relative  and  absolute 
orientation  are  found  simultaneously  with  a  so-called  bundle 
adjustment.  In  our  system  this  approach  is  the  one  used  for 
radar  images.  The  basic  equations  for  the  radar  stereo 
adjustment  consist  of  two  types: 

(a)  The  squint  angle  condition  (see  E  Q  .  2.10): 

£•(£-£)  -  sin:-  s  •  £-_s  =  0 

(b)  The  range  condition  (see  EQ  .  2.9): 

r  s  "  P-1'  =  0 

for  slant  range  presentations,  or 

V"r  :  +  -  p  -  s  =0 

g  r  - 

for  ground  range  presentations. 

Each  measured  point  gives  rise  to  two  equations  of  type 
(4.1)  and  (4.2).  In  addition  each  pair  of  homologue  image 
points  produces  one  additional  condition.  This  is  called 
the  cocircularity  condition,  which  defines  the  location  of 


(4.2a) 


(4.2b) 
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an  imaged  point  in  object  space  by  the  intersection  of  two 
circles,  each  defined  by  two  equations  of  type  (4.1)  and 
(4.2).  The  equations  are  r.or.-lii ear.  Linearized  forms  are 
used  in  an  iterative  solution. 

Observations  for  the  computations  are  sensor  position 
measurements,  image  coordinate  measurements  of  ground 
control  points  in  one  or  in  both  images  and  pairs  of 
homolog ue  image  coordinate  measurements  of  stereo  model 
orientation  points.  Unknowns  to  be  determined  within  the 
adjustment  are: 

-  the  coefficients  of  the  orbit  time  polynomials; 

-  the  radar  imaging  parameters  (inner  orienta¬ 
tion)  to  convert  radar  image  coordinates  to 
range  and  time; 

-  the  parameters  of  a  correction  polynomial  to 
describe  radar  image  Reformations; 

-  a  value  for  the  squint  angle. 

The  solution  for  the  unknowns  i r.  each  iteration  of  the 
non-linear  equation  system  is  obtained  in  the  instruments 
host  processor  in  a  least  squares  adjustment  using  t'ne 
method  of  conjugated  gradients  (Schwarz,  19  7  0)  .  This 
solution  method  is  iterative  and  may  have  sere  advantages 
for  larger  equation  systems  when  a  small  computer  must  be 
used  and  both  computing  times  and  storage  requirements  are 
limited.  Furthermore  there  is  evidence  that  the  method  is 
more  robust  for  this  special  problem  than  a  strict  Gauss 
solution  of  normal  equations. 

Depending  on  the  quality  of  the  approximate  values  the 
adjustment  process  must  be  repeated  iteratively  until 
results  satisfy  specified  termination  criteria. 


i 
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4.4  Example  J  f  Radar  3_tereo  Pl^oJ^t^ng  ££  SIR - A  Satellite 

Data 

A  Space  Shuttle  Imaging  Radar  SIR-A  stereo  image  pair 
of  the  Greek  islands  Cephalonia  and  Ithaka  was  used  to 
illustrate  the  procedure  of  extracting  a  digital  height 
model  (Fig.  4.1).  A  more  detailed  study  of  the 
stereoscopic  computations  from  such  data  is  presented  by 
Kobrick  et  a  1 .  (submitted).  Viewing  limitations  were 
discussed  by  Domik  et  al.(1983).  After  model  set-up  the 
topographic  elevations  were  digitized  in  the  object 
coordinate  system  using  profiles  and  ridge  lines;  the 
measurements  were  entered  into  the  terrain  model  system  G  T  M  . 
The  GTM  program  enables  one  to  generate  contour  line  plots 
at  chosen  interval,  to  obtain  axonometric  views  of  the 
terrain  data  and  various  other  products. 

Figure  4.2  shows  an  axonometric  view  of  the  digital 
elevation  model  (DEM)  of  the  islands  extracted  from  a  map 
1  :  200  000,  Figure  4.3  presents  the  DEM  created  on  the  Kern 
DSR-1  analytical  stereoplotter.  A  comparison  of  the  two 
OEM's  is  possible  with  the  GTM-system  and  reveals  that  the 
root  mean  square  differences  over  all  DEM  points  of  the 
islands  amount  to  +  98  meters. 


37 


r 


FIGURE  4.1:  Space  Shuttle  Imaging  Radar  SIR-A 

stereo  radar  image  pair  of  the  Greek 
islands  Cephalonia  and  Ithaka,  taken 
at  34°  angle  between  flight  lines 
at  45°  look  angle  off-nadir. 
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FIGURE  4.3:  Digital  elevation  model  of  the  greek 
islands  Cephalonia  and  Ithaka  in 
axonoraetric  view.  Data  digitized  at 
KERN  DSR-1  stereoplotter  using 
program  SMART. 
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4.5  Software  Implementation 


The  software  package  is  currently  installed  at  a  Kern 
DSR-1  analytical  plotter  under  the  name  SMART  (Stereo 
Mapping  with  Radar  Techniques). 

The  DSR-1  standard  configuration  consists  of  three 
independent  digital  processing  units,  which  are  denoted  as 
PI,  P  2  and  P3  (Chapuis,  1  9  8  0  )  .  Processor  PI  is  the  host 
computer.  It  serves  for  development  and  execution  of 
application  programs  and  performs  data  transfer  between  the 
other  processors.  The  computation  of  set-up  parameters  for 
the  radar  stereo  model  is  done  in  PI.  Processor  P2  receives 
the  radar  image  orientation  parameters  from  PI  and  uses  them 
to  convert  in  real  time  model  coordinates  X,Y,Z  to  radar 
image  coordinates  r ' , t ' , r " , t "  and  onto  DSR-1  plate 
coordinates  x',y'  and  x",y".  P2  also  transfers  image,  model 
or  object  coordinates  to  processor  PI  if  requested. 
Processor  P3  is  an  added  convenience  for  the  operator  and 
controls  an  operator  control  panel  which  may  be  used  for 
communication  between  operator  and  DSR-1. 

The  Kern  DSR-1  therefore  uses  three  computers  as 
opposed  to  other  plotters  which  may  employ  only  one.  What 
seems  to  be  an  added  complexity  in  fact  is  an  advantage: 
the  real-time  operations  are  performed  by  a  seperate 
processor  P2  with  a  clearly  defined  software  interface  for  a 
user  working  with  the  host  PI.  The  user  will  therefore  work 
in  a  familiar,  general  purpose  multi-user  environment  with 
an  unmodified  operating  system.  The  critical  real  time 
operations  are  singled  out. 


i 


5 . J  RADAR  IMAGE  SIMULATION 

Approaches  to  image  simulation  can  be  categorised 
according  to  the  domain  of  the  transformation  function. 
Object  space  algorithms  map  each  object  point  (DEM  point) 
into  a  two  dimensional  field  (image)  whereas  image  space 
algorithms  start  from  equidistant  image  coordinates  (i,j)  of 
the  output  image  and  assign  gray  value  accordingly  to  it. 
Object  space  algorithms  are  usually  faster  in  their 
application  than  image  space  algorithms.  But  there  is  a 
need  for  interpolation  to  convert  the  unequally  spaced  image 
addresses  to  equidistant  ones. 

The  development  of  the  radar  simulation  was  part  of 
another  project  (see  Domik,  Leberl  and  Raggam,  1983).  It  is 
of  great  relevance  in  the  current  context  and  is  discussed 
in  the  following  for  ease  of  reference. 

An  image  space  algorithm  was  chosen  to  be  implemented 
because  it  provides  DEM  addresses  at  exact  image  locations. 
Fig.  5.1  and  5.2  show  the  basic  idea  of  object  vs.  image 
space  algorithms. 

Gray  values  are  attached  to  the  image  coordinates 
dependent  on  the  corresponding  incidence  angle  on  the 
ground . 


5 • 1  The  leone t rio  Imaging  Model 

The  image  coordinates  (i,j)  serve  as  start  values  for 
the  simulation  program.  They  are  used  to  extract  imaging 
time  t  and  slant  range  r  or  ground  range  g  by  Eq.'s  5.1. 


t  = 

(i-1 )*f . 

l 

+ 

C* 

O 

EQ.  (5.  la) 

r  = 

(j-n*fj 

+  3d 

EQ. (5.  lb) 

g  = 

(j-n*fj 

+  sd 

r  = 

SQRT(H2 

♦  s2) 

EQ. (5.1c) 

where 
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The  sensor  position  and  velocity  vector, 
derived  from  the  imaging  time,  define  a  cartesian 
coordinate  system:  The  sensor  position  is  the 

origin  of  this  "sensor  coordinate  system",  the  axes 
are  defined  by  the  velocity  vector  and  the  sensor 
position  as  expressed  in  a  geocentric  coordinate 
system  (Eq.  2.8).  The  set-up  of  this  system  allows 
a  clear  presentation  of  the  geometry  created  by  any 
input  coordinate  (i,j). 

The  origin  of  the  sensor  coordinate  system  is 
the  center  of  a  sphere  with  radius  r  (range  sphere). 
In  general,  a  squint  angle  t  i  0  is  assumed  for 
imaging.  Thus  a  cone  is  defined  with  the  peak  in 
the  origin  and  its  axis  along  the  velocity  vector  of 
the  air-  or  spacecraft.  All  generating  lines  form 
an  angle  of  (90°  -  -  )  with  the  axis.  If  the 
squint  angle  t  is  equal  to  0  imaging  is  in  a  plane 
perpendicular  to  the  velocity  vector.  The 

intersection  of  sphere  and  cone,  or,  sphere  and 
plane,  respectively,  defines  a  circle  (intersection 
circle).  The  remaining  task  is  now  to  intersect  the 
circle  with  the  DEM. 

At  least  one  intersection  point  but  possibly 
several  ones  will  be  found.  The  occurrence  of  more 
than  one  solution  is  called  "lay-over":  different 

terrain  targets  were  illuminated  simultaneously. 
The  address  of  each  intersection  point  as  well  as 
corresponding  look  angle,  surface  inclinations  and 
additional  informations  are  stored  for  further 
processing. 


5.2  Radiometric  Model 


The  definition  of  an  image  gray  value  in  a 
given  pixel  location  is  denoted  as  "radiometric 
model".  One  has  to  define  areas  of  shadow  (black), 
confuse  reflection,  multiple  reflection  in  lay-overs 
etc. 


Black  image  areas  are  mostly  shadows  or  water 
bodies.  In  order  to  find  those  pixels  that  indeed 
are  shadow  areas  one  has  to  sort  the  look  angles 
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going  from  near  range  to  far  range  along  an  imaged 
object  line.  Shadow  areas  begin  where  an  increasing 
range  r  is  associated  with  a  reduced  look-angle  off 
r.adir. 

In  areas  other  than  shadows  the  image  gray 
value  is  a  function  of  the  incidence  angle  for  the 
incoming  radiation.  This  is  an  angle  between  the 
line  of  sight  antenna-object  and  the  local  normal 
onto  the  surface. 

Once  this  angle  is  computed  for  a  pixel  one  has 
the  choice  of  selecting  a  back  scatter  curve 
according  to  the  thematic  charactaristics  of  the 
material  and  of  the  reflection  angle  t  .  Standard 
back scatter  function  are  those  of  Hagfors  (1964), 
Muhleman  (  1  964  )  or  a  c os i n e- f unc t  i  o n  . 

Simultaneous  illumination  of  various  terrain 
parts  leads  to  image  lay-overs:  the  projection 
circle  intersects  more  than  one  point  on  the  ground. 
? h  •  gray  value  is  the  sum  of  the  individual 
reflection  values. 


5 . 3  Data  Structure 

The  internal  data  structures  are  constructed  to 
make  available  a  fast  access  to  the  address  file 
elements  in  relation  to  the  image  elements.  Each 
3 z i m u t h  line  in  the  image  represents  a  logical  block 
and  corresponds  to  a  unique  logical  block  of  the 
address  file.  Memory  allocation  for  each  logical 
block  will  be  a  minimum  cf  r.  *  e  bytes,  where  r.  is  the 
number  of  range  columns  and  e  the  number  of  bytes 
used  for  storage  of  one  gray  value.  The  address 
file  allocates  (at  minimum  again)  h  *  f  bytes,  with  h 
being  the  number  of  addresses  belonging  to  one 
azimuth  line  and  f  the  number  of  bytes  used  for  one 
address  definition.  The  relation  is  realized  by  an 
additional  a c c e s s -ma t r i x  (pointer  file),  which 
allows  a  fast  associatior  between  elements  of  the 
two  different  files.  The  additional  storage  for  the 
relation  matrix  will  be  h  *  2  bytes  per  logical  block. 
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5.4  Software  Implementation 

The  product  of  the  software  implementation 
perform  the  simulation  is  named  3 1 M  R 1 3  A  ( 3  I  M u 1 a  t  i  :  n 
of  Radar  images  using  an  loage  Space  Algorithm  <  a r.  : 
is  currently  running  on  a  7  AX-11/750. 

Emphasis  was  put  on  modularity  of  the  system 
(see  description  of  system  modules,  Fig.  5.3).  The 
program  can  be  easily  changed,  e . g .  additional 
thematic  input  data  sets  could  allow  for  use  :  f 
different  backseatter  curves  or  use  of  different 
look-up  tables  within  one  DEM  area. 


6  .  0 


3THER  AUXILIARY  SYSTEM  ELEMENTS 


o . 1  Registration 

The  image  registration  is  based  on  a  program 
RECTI?  by  Diarra  (1982);  it  is  a  modular  software 
system  for  geometric  correction  of  an  image.  Input 
is  a  set  of  homologue  pairs  of  points  whose 
coordinates  are  known  in  the  real  image  and  in  the 
simulation.  These  are  used  to  compute  correction 
polynomials  to  describe  the  trend  of  the  image 
deformations.  Residuals  at  control  points  serve 
then  to  compute  refined  corrections  with  some  higher 
order  i.nterpoiative  method.  Work  is  thus  via 
"warping"  or  "rubber-sheeting".  Polynomials  are 
limited  to  degree  3  or  less,  so  that  they  are  both 
stable  and  overdetermined.  Residuals  are  extracted 
at  each  anchor  point,  and  serve  to  mod  1 fy  the 
original  polynomials. 

These  correction  polynomials  are  used  for 
resampling  based  on  a  grid  of  anchor  points;  each 
set  of  object  coordinates  obtains  image-positions. 
Various  methods  may  then  be  applied  to  determine 
gray  values  at  those  positions:  interpolation 
followed  by  nearest  neighbor  assignment,  bilinear 
interpolation,  and  bicubic  interpolation,  are  the 
current  options.  Control  points  will  usually  be 
defined  by  hand,  or  by  digital  correction. 

The  registration  results  in  a  rectified  image; 
in  the  case  of  data  actually  studied  remaining 
distortions  are  of  an  average  of  3  pixels  or  less, 
assuming  an  appropriate  distribution  of  control 
points  . 

This  results  in  a  real  radar  image  superimposed 
onto  the  simulated  data. 
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capability  of  the  corresponding  targets.  With  fixed 
wavelength  and  polarisation  the  brightness  varies 
only  with  incidence  angle  and  reflection  properties. 


The  location  of  the  pixel 
defined  by  the  time  delay  of  the 
from  the  target.  For  a  perfectly 
time  delay  increases  proportional 
nadir  and  target.  Irregularities 
be  imaged  affect  the  time 
distortions  in  the  image. 


in  the  image  is 
radar  pulse  to  and 
flat  ground  the 
to  the  distance  of 
on  the  surface  to 
delay  and  create 


An  extreme  effect  of  relief  displacement  occurs 
with  lay-over;  the  top  of  a  target  to  be  imaged  is 
closer  to  the  radar  sensor  than  the  bottom  and  is 
therefore  recorded  sooner.  However,  lay-over  is 
also  defined  as  the  simultaneous  mapping  of  two 
different  ground  targets  into  one  resolution  cell. 
This  is  the  case  when  the  slant  range  of  two  objects 
is  the  same.  Lay-over  effects  the  brightness  as  the 
return  power  from  two  or  more  targets  is  added  up. 
The  simulation  process  can  correctly  consider  the 
lay-over  case  by  attaching  more  than  one  position 
address  to  an  image  gray  value. 

After  the  simulation  and  registration  each 
digital  D  N  number  of  the  corrected  radar  image  is 
connected  to  one  or  more  DEM  addresses.  The 
procedure  GEO  R  EC  (GEOmetric  RECtification)  performs 
the  final  step  of  the  radar  image  rectification 
making  use  of  the  data  structures  described  in 
chapter  5.3. 


Input  to  this  last  step  are 

-  grey  values  of  the  real  image  registered 
to  the  simulation 
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address  file  (containing  DEM  addresses) 


-  relation  between  grey  values  and  address 
file 

According  to  the  data  structure  description  the 
logical  blocks  of  grey  value  and  address  file 
correspond  sequentially  to  each  other.  This  allows 
to  process  each  azimuth  line  ( corresponding  to  one 
logical  block)  independently.  Each  gray  value  is 
now  placed  at  the  DEM  addresses  calculated  during 
the  simulation.  Thus  the  geometric  radar 
distortions  are  eliminated  again  and  the  resulting 
ortho  image  displays  distances  instead  of  (slant  or 
ground )  ranges . 

Usually  the  DEM-addresses  are  not  exact  grid 
positions,  so  that  some  kind  of  interpolation  has  to 
take  place. 
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7.0  EXAMPLES  OF  RADAR  IMAGE  RECTIFICATION 


7 . 1  Oetztal-Data 

A  radar  image  was  obtained  during  the  European 
S A R - 5 3 0  campaign  of  the  Oetztal  in  Tyrol  (Rott, 
1983).  The  sensor  for  this  mission  was  carried 
aboard  a  Convair-580  aircraft  with  a  dual  band,  dual 
polarisation  SAR  which  could  be  operated  in  X  and  C 
band  or  X  and  L  band.  Data  were  recorded  by  using 
either  an  optical  or  a  digital  recorder: 

The  aim  of  the  Tyrol-part  of  the  SAR -530 
experiment  were  investigations  in  the  applications 
of  radar  for  mapping  snow  and  glaciers  in 
mountainous  regions.  The  sensor  was  operated  in  X 
and  C  band  (HH  polarisation),  covering  a  ground 
swath  of  13  km  length.  Fig.  7.1  is  a  presentation 
of  the  SAR  image,  X  band,  recorded  digitally  and 
processed  at  the  German  Aerospace  Research 
Establishment  ( D F V L R  )  .  The  presentation  is  in 
ground  range;  the  original  resolution  of  3  m  was 
changed  to  a  coarser  resolution  of  2d  m  for 
presentation  in  512  x  512  pixels. 

A  DEM  was  derived  from  a  map  1  :  50  000  by 

scanning  the  contour  and  drain /ridge  lines.  The 
imaged  part  is  only  a  subarea  of  the  DEM,  which 
extends  from  1  0  0  4  8  '  to  10°59'  East  longitude  and 
from  4  6  0  4  5  '  to  4  6  0  5  2  ’  North  latitude.  The 
mountain  spreading  from  the  south-western  towards 
the  nor th-eastern  corner  can  be  recognized  in  the 
radar  image  (upper  half).  The  illuminated  digital 
height  model  of  the  area  to  be  investigated  is 
presented  in  Fig.  7.2. 
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7 . 2  Processing  Of  The  Oe  tz  ta  1  Data 


The  simulation  was  carried  out  by  the  simulator 
according  to  chapter  5.  Since  flight  navigation  was 
eroneous  during  the  SAR-580  flight,  one  has  to  rely 
on  control  points  to  approximately  derive  the  flight 
and  imaging  parameters.  The  maximal  height  in  this 
area  is  3615  m,  whereas  the  minimal  height  is  85  m. 
Hagfors  law  was  used  to  model  the  backscatter 
depending  on  the  incidence  angle.  No  additional 
input  was  used  to  distinguish  between  different 
thematic  features  in  the  image. 

The  simulation  results  in  a  gray  value  image 
with  an  additional  address  file  and  relation  matrix. 
The  gray  values  assigned  by  the  simulator  are  merely 
functions  of  the  incidence  angle  and  do  not  consider 
any  properties  of  the  imaged  targets  (Fig.  7-3). 

The  geometry  created  by  the  simulation  matches 
the  real  image  closely.  The  radar  distortions  due 
to  radar  properities  like  forshortening  and 
shadowing  are  correctly  implemented  by  the 
simulation . 

However,  sensor  dependent  errors  will  be 
present  due  to  the  inaccurate  flight  and  imaging 
records.  Thus  there  is  a  need  for  geometric  image 
correction.  The  registration  described  in  chapter  6 
creates  a  modified  radar  image  to  fit  with  the 
simulated  one.  The  simulated  image  serves  to 
support  control  point  definition. 

The  corrections  applied  to  the  real  radar  image 
are  shown  in  Fig.  7 . 4 .  Corresponding  features  in 
both  real  and  simulated  image  were  used  as  control 
points  in  the  registration.  Their  residual  vectors 
are  plotted  into  the  deformation  grid  with  a  10  x  10 
pixel  size.  The  DN  values  of  the  new  SAR  image 
(Fig.  7.5)  may  now  replace  the  synthetic  gray 
values  and  correspond  to  the  additional  output  files 
of  the  simulation.  Thus  the  connection  between 
radar  image  and  DEM  addresses  is  made  and  the 
rectification  can  be  carried  out  by  the  appropriate 
procedure. 
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The  radar  ortho  image  is  shown  in  Fig.  7.6. 


IGURE  7.3s  Simulated  radar  image:  approximated 
flight  path  and  elevation  model 
serve  as  input.  Ground  range  presen 
tation  and  diffuse  scattering  are 
assumed  . 
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The  simulation  displays  the  gray  values  as  a 
function  of  the  incidence  angle,  whereas  the  gray 
values  in  the  real  image  vary  with  the  change  of 
incidence  angle  and  thematic  characteristics. 
Substracting  the  simulated  gray  values  from  the  real 
data  leaves  a  radiometric  corrected  image  (Fig. 
7.7):  Those  parts  in  the  image,  which  are  mainly 
influenced  by  the  surface  inclination  are  displayed 
in  middle  gray  tones,  parts  where  the  thematic 
content  enhances  the  return  power  are  to  be 
distinguished  by  a  lighter  gray  to  white.  Dark  gray 
tones  indicate  the  strong  reflection  in  the 
simulation  as  compared  to  the  real  image. 
Especially  near  the  edges  of  the  image  the  dark 
tones  are  largely  insignificant. 


FIGURE  7.7:  R ad i ome t r i ca 1 1 y  corrected  radar  image 
Influence  of  topography  in  radar 
brightness  is  removed. 


» 

7 . 3  Cephalonia  Data 


On  Nov.  12,  1981,  the  SIR-A  (Shuttle  Imaging 

Radar-A)  was  launched  aboard  NASA's  space  shuttle. 
SIR-A  is  a  left-side-looking  SAR  operated  at  L-band 
frequency  with  an  elevation  angle  of  approximately 
50°  off  nadir  to  the  center  of  the  imaging  swath. 
SIR-A  data  were  optically  recorded. 

Two  images  obtained  in  Europe  from  Data  Take 
32-33  and  Data  Take  37 A,  resp .  ,  were  of  special 
interest  for  stereo  investigations  (see  Fig.  4.1). 
Though  SIR-A  was  not  specially  designed  for 
stereo-investigations,  several  data  takes  from 
crossing  orbits  proved  to  deliver  fairly  well 
stereo.  The  intersection  angle  of  the  imaged  swaths 
was  approximately  34°. 

For  our  investigations  the  images  presented  on 
film  were  scanned  by  an  Optronics  device  to  receive 
the  data  in  digital  form.  The  digital  data  was 
processed  to  receive  a  resolution  of  approximately 
100m. 


The  analog  images  were  used  to  create  a 
stereo-derived  elevation  model  as  explained  in 
chapter  4.  At  the  same  time  an  elevation  model  was 
genereated  by  digitizing  the  contour-lines  at  a  200 
m  interval  of  a  1  :  200,000  map  in  UTM  projection. 
Axonometric  presentations  of  both  models  are 
presented  in  Fig.  4.1  . 


7  .  4  Processing  0 f  Cej>hal_onia  Da _ta 


The  process  is  analog  to  the  Oetztal  data 
processing,  using  the  modules  of  the  radar  image 
processing  system  described  in  chapters  4  to  6. 
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Flight  and  imaging  information  to  be  entered 
into  the  simulation  were  gained  from  (Cimino  and 
Elachi,  1982):  Latitude  and  longitude  for  the 
center  of  the  swath  as  well  as  flight  altitude  were 
recorded  only  for  each  minute,  so  the  flight  path 
was  approximated  by  a  straight  line  at  the  height  of 
about  264  km  (slightly  descending).  The  STC 
position  and  the  interpulse  period  were  used  to 
calculate  the  sweep  delay  to  the  near  edge  of  the 
swath.  The  image  presentation  was  in  slant  range. 
Through  the  simulator  a  synthetic  radar  image  was 
created  from  the  stereo-derived  DEM  together  with  an 
additional  address  file  (Fig.  7.3).  The 
registration  resulted  in  Fig  7.9  bv  removing  the 
differences  (see  registration  grid  in  Fig.  7.10) 
between  real  and  simulated  image.  The  final  ortho 
image  (Fig.  7.11)  was  created  by  resampling  the 
radar  image  grey  values  at  their  corresponding  DEM 
addresses . 


'IGURE  7.8:  Simulated  radar  image  using  imaging  and 
sensor  parameters  from  the  5IE-A  record  i 
and  stereo-derived  elevation  model. 


CONCLUSION  AMD  OUTLOOK 


3  .  0 


We  have  described  an  elaborate  software 
for  the  rectification  ot  digital  radar 
Rectification  is  based  on  digital  elevation 
from  maps  or  from  radar  stereo  models,  on 
points  and  on  given  data  of  the  sensor  and 
path. 
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The  core  of  the  procedure  is  an  image 
simulation  to  support  the  identification  of 
homologous  features  in  both  the  image  and  the  DEM. 
Rectification  can  be  geometric  and/or  radiometric. 


The  entire  system  is  experimental:  the  size  of 
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